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O !■ Abstract. Wc consider a simple decision model in which a set of agents randomly 

q | choose one of two competing shops selling the same perishable products (typically 

O ' food). The satisfaction of agents with respect to a given store is related to the freshness 

o 



of the previously bought products. Agents select with a higher probability the store 
they are most satisfied with. Studying the model from a statistical physics perspective, 
both through numerical simulations and mean-field analytical methods, we find a rich 
behaviour with continuous and discontinuous phase transitions between a symmetric 
i-C ' phase where both stores maintain the same level of activity, and a phase with broken 

symmetry where one of the two shops attracts more customers than the other. 



1. Introduction 



> 
oo 

bi- 
phase transitions and symmetry breaking are standard phenomena occuring in physical 
O systems due to the presence of interactions between the involved microscopic entities 

(most often atoms or molecules) p]. A prominent example is the ferromagnetic Ising 
model [2], in which microscopic magnetic moments, schematically represented as spin 
variables Sj = ±1, get aligned at low temperature due to the interaction energy that 
couples neighbouring spins. Outside physics, phase transitions are also often found in 
complex systems composed of many interacting "agents" (virtual entities representing 
in a schematic way human beings, animals or bacteria for instance) [3j HI [5]. In these 
cases, interaction rules do not derive from an energy function as in physical systems, 
but still correspond to direct interactions between agents. 

There however exist interesting situations in which agents have no direct 
interactions between them, but rather interact with an external system, having its own 
dynamics, and producing a feedback on the agents. A simple example of such a situation 
is the case of an assembly of customers who go for shopping on a regular time basis -say, 
every week- and need to choose between several equivalent stores. In the absence of 
bias between shops, that is if customers choose any of the shops with equal probability, 
all stores should maintain the same activity level. If however, a fluctuation of activity 
occurs, one shop may become slightly less attractive because, for instance, the sold 
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products get older (which may have significant consequences when selling, e.g., fresh 
food). One can intuitively expect that this feedback loop may, in some cases, result in 
an amplification mechanism leading to a symmetry breaking between shops. 

Modelling the competition between stores is a difficult issue, which raised interest in 
the economic literature from different perspectives [El El El El \10\ as well as, marginally, 
in the physics literature [ITj . Such models are often grounded in discrete choice theory 
[T2~l [T3j . Despite the central role played by prices in economic modelling, empirical 
studies show that pricing mechanisms for fresh products like fruits and vegetables are 
complex, and that other criteria related to the quality of the products also play an 
important role in customers decisions [H]. The freshness of products thus turns out to 
be an important ingredient to take into account when modelling shops selling perishable 



In the present paper, we investigate from a statistical physics viewpoint a simple 
model in which the ageing of products, and the corresponding decrease of utility, are the 
central ingredients. At variance with economic models, we neglect all financial aspects 
(like the price of products and the fixed costs that the shopkeeper has to pay), in 
order to single out a unique dynamical feedback mechanism (the freshness of products) 
between shops and customers, thus leading to a minimal dynamical model of shop 
competition and decision making. Focusing on the symmetry breaking phenomenon, we 
find both through numerical simulations and analytical mean-field approaches that a 
phase transition takes place between a symmetric state in which the two stores maintain 
the same level of activity, and a symmetry broken state where a majority of customers 
choose a given store, in which they can find fresher products. 

2. Model of competing stores 

2.1. Definition 

The model is defined as follows. A set of N a agents randomly visit, on a regularly time 
basis, one of two competing stores selling the same products. Each store has the same 
number N p of products, which are all identical except for their age. The age r defines 
the freshness h(r) of the corresponding product, where h(r) is a decreasing function, 
taking values between and 1, of the variable r > 0. A typical functional form for h{r) 
is 



where T\ is the characteristic ageing time of products. 

At discrete time steps tk = k5t, with 5t = T Q /N a , an agent is chosen at random 
with uniform probability. This agent i then randomly selects one of the two stores with 
a probability related to its satisfaction SV,- with respect to each store j = 1,2. To be 
more specific, the probability Pi that agent i chooses the first store is given by 



goods [ZIIHICEE]. 




(1) 



1 



(2) 



Pi = 



1 + e-ASi/r 
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Figure 1. Illustration of the dynamical behaviour of the model, in cases when 
symmetry breaking occurs (left) or does not occur (right); the mean freshness (hj) 
of products in each shop j = 1, 2 is plotted as a function of time. The inset on the left 
plot shows a short time- window corresponding to the transient state of the simulation 
starting from a non-equilibrated initial conditions. Left panel: t\ = 16.7; right panel: 
Ti = 40. The other parameters are the same on both panels: N a = N p — 500, T = 0.1, 

T = l. 

where ASi = Sn — Si2 is the satisfaction difference of agent i with respect to the two 
shops, and T is a parameter playing the role of a temperature, taking into account the 
influence of other factors that are not described explicitly. The second store is chosen 
with probability 1 — pj. This dynamical rule preserves the symmetry between the two 
shops. Note that Eq. (j2J) is known as the logit rule in the socio-economic literature [32] , 
and that it is also formally similar to the Glauber transition probability used in physics, 
provided —Si can be interpreted as an energy (which is not necessarily the case in social 
models, see e.g. [16]). 

Once agent i has chosen a store j = 1,2, it visits this store and buys one product 
chosen at random (we recall that only one type of product is sold, only the age r differs 
from one to the other). The sold product is immediately replaced by a fresh one, with 
age t' = 0. The agent a posteriori assesses the freshness h{r) of the bought product, 
and updates its satisfaction with respect to the store j according to 

S'( 7 . = a^ + (l-a)/i(r), (3) 

where < a < 1 is a parameter characterizing the agents' memory (all agents have the 
same value of a). Then the process is iterated and a new agent is chosen randomly. 

2.2. Qualitative behaviour of the model 

Simulating numerically the above model, we observe for low enough temperature a 
symmetry breaking between the two stores (see the left panel of Fig. [TJ): one of the 
stores reaches a high level of commercial activity, resulting in a high freshness of the 
sold products, while the second one only maintains a low activity and a low freshness 
level. Interestingly, symmetry can be restored by varying the parameter n, which 
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controls the ageing rate of products, as seen on the right panel of Fig. [TJ both shops 
then keep the same level of freshness. Less surprisingly, symmetry is also restored at 
high temperature (not shown). 

In the following sections, we quantitatively study this symmetry breaking phase 
transition, through both analytical and numerical methods. Section [3] is devoted to a 
mean-field approach that allows us to obtain the phase diagram of the model. We find 
in particular a transition line in the plane defined by the temperature T and a control 
parameter R related to the characteristic ageing time T\. Section H] is then devoted to 
the comparison between mean-field results and numerical simulations. Finally, Section [5] 
briefly discusses an analogy with the Ising model, one of the most standard model 
exhibiting symmetry breaking. 



3. Mean- field approach 

3.1. Framework and basic dynamical equations 

In this section, we present an analytical approach of the above model, in order to 
characterize quantitatively the phase transition between symmetric and non-symmetric 
states of the two stores. Our analytical framework relies on the simplifying assumption 
that the fluctuations of the agents' satisfaction can be neglected, so that all agents are 
considered to have a satisfaction difference AS* equal to the average value AS. In socio- 
economic terms, one would say that the set of agents is replaced by a "representative 



agent" [12] . To lighten the notations, we introduce a = AS. From Eq. ([3]), we get, 
assuming that the time spent between two visits to a shop is exactly tq (while this is 
only true on average): 

a{t + r ) = aa(t) + (1 - a) (h^t) - h 2 (t)) , (4) 

where hj(t) is the average freshness of products in the store j at time t. This equation 

can be reformulated as 

1 1 — n> 

- (a(t + r ) - ait)) = (-a(t) + h^t) - h 2 (t)) . (5) 

r TO 

Hence <r(t) varies significantly over a time scale of the order of t /(1 — a). In the limit 
of strong memory, < 1 — a <C 1, r /(l — a) ^> r so that the variations of a are 
very small over a time To Hence <r(t) can be considered as a continuous function of the 
variable t, and the l.h.s. of Eq. §5§ can be replaced by a time derivative. We thus end 
up with 

^= 7 -£,(*)), (6) 

with 7 = (1 — a)/ro. It is convenient to rather express the dynamics in terms of the 
fraction p of customers choosing the first store, instead of keeping the variable a. Using 
Eq. ([2]) and replacing ASi by the average value a, we get the relation 



Phase transition in a dynamical decision model 5 
which can be plugged into Eq. fl§]), leading to 

(8) 



tt=f P{1 - p) 



V 



Tin I ^— ^ ) +h 1 (t)-h 2 (t) 



We now need to compute the average freshness hj as a function of p. To this aim, we 
first determine the age distribution of products. 

3.2. Age distribution of products 

Let us denote as \j(t) dt the average fraction of products sold in store j in the time 
interval [t,t + dt], and <f)j(r,t) the probability distribution of the age r of products in 
store j at time t, normalized according to J °° (f>j(r,t) dr = 1. This distribution evolves 
in time according to the following relations: 

0j(r + dt,t + dt) = (f>j(r, t) - Xj(t) dt ^(r, t) for r > 0, 
4>j(0,t) dt = Xj(t) dt. 

The second equation in ([9]) encodes the fact that sold products, selected at random, are 
immediately replaced by new ones with age r = 0. Expanding the first equation in ([9]) 
to linear order in dt, one finds 



r (r,t) + ^(T,t) = -X ] (t) ( f> J (T,t). (10) 

Note that integrating this equation over r > yields 0,(0, t) = Xj(t), consistently with 
the second equation in Qj. The stationary solution of Eq. (flOj) is exponential, namely 

#(t) = V"V. (11) 

Using the exponential ageing function given in Eq. (CQ), the corresponding average 
freshness is thus 

h 3 = dT<pf(r) = (l + J-) . (12) 

In addition, the fraction Xj dt of products sold in the time interval [t, t + dt] is obtained 
as the ratio of the number N a pj dt/ro of customers visiting shop j in this interval, and 
of the number N p of products, yielding 

x > - W» (13) 

with pi = p and p 2 = 1 —p- We thus end up, as far as the average freshness is concerned, 
with 

where we have defined the control parameter 
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Figure 2. Illustration of the behaviour of the free energy /T,i?(m), for four different 
temperatures (from top to bottom: T = 0.2, 0.185, 0.16, and 0.12) and the same value 
R = 0.2, showing the onset of the two symmetric secondary minima. 



The dimensionless ratio R compares N p , the total number of products available in a shop, 
to ^ T i, the typical number of agents visiting a shop during the characteristic ageing 
time of the products (assuming a symmetric state). A small value of R thus corresponds 
to a situation where products are often renewed, while for large R products are rarely 
selected and become old on average. In the following, we use as an approximation the 
above steady-state distribution of the age of products, assuming that the satisfaction 
of agents evolves on time scales much longer than the turn-over time of products. This 
assumption has, however, no consequence on the determination of the fixed points of 
the dynamics. 



3.3. Determination of the phase transition 

Gathering results, we find for the time evolution of p the following equation 



dp 
It 



7 



p(l-p) 



Tin 



V 



V 



R 

2p 



-i 



R 



2(1 -p) 



-i 



.(16) 



Introducing an order parameter m through p — (1 + m)/2, we can formally rewrite 
Eq. flUD as 



dm ^ 2 

= ? (l-m) 9T , R (m), 



with grp^m) given by 

9T,R( m ) : 



Tin 



1 — m 
1 + m 



R 



m 



-i 



R 



1 — m 



-i 



;i7) 



.(18) 



The reason for excluding the factor 7(1 — m 2 )/T from the definition of grfi^) will 
be discussed in Sec. El In order to study possible phase transitions, it is convenient to 




Figure 3. Mean-field phase diagram, with four different regions. A: m = is the 
only minimum of fT.Ft.{ m )- B: m = is the global minimum, but secondary minima 
have appeared. C: m = is only a local minimum, the two symmetric minima are the 
most stable. D: m = is a local maximum, the two symmetric minima are the only 
stable states. The transition lines between those regions are respectively defined by 
T(R) = T s between A and B, T<j between B and C, and T c between C and D. T c also 
separates A and D for R > R* ps 0.732, when the transition is continuous. The lower 
panels show the typical shape of fT,ii( m ) for m > {fT,R( m ) is a symmetric function) 
in each region. 



define an effective free energy (or potential function) /r,i?(w) through 



so that 



dm 
~dt ~ 



9T,R{ m 



dm', 



df: 



T.R 



dm 



(19) 



(20) 



It is easily seen that m = 0, corresponding to the symmetric state, is a steady-state 
solution of Eq. ( fTTl) . To test the local stability of this solution, we expand the free 
energy fT,R( m ) ^ or sma ^ m i to order m 2 , leading to 



fT,R( m ) 



R 



m 



0{m 4 



(21) 



(1 + R) 2 \ 

The solution m = is linearly stable if the bracket on the r.h.s. is positive. The critical 
temperature T c is reached when the bracket vanishes, yielding 

R 



^ C ~{l + Rf (22) 
When this transition coincides with the onset of two symmetric minima for T < T c , 
the transition is continuous. For low values of R, however, two symmetric secondary 
minima appear already for a value T s > T c , when m = is still the global minimum as 
illustrated on Fig. [2J This case corresponds to a discontinuous transition, which should 
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Figure 4. Order parameter m, plotted either as a function of R or as a function of 
temperature T. Lines correspond to the mean-field predictions, while dots with error 
bars are the results from numerical simulations. On all four figures, a phase transition, 
either continuous or discontinuous, is observed between a symmetric phase (both stores 
have the same activity) and a phase with broken symmetry (one of the stores has a 
higher activity level). 

be reached (in analogy to the equilibrium theory of phase transitions) at a temperature 
Trf such that all three minima have equal values fT,ii{ m ) = 0. 

The resulting mean-field phase diagram is presented on Fig. |3] In region A, a single 
steady-state solution (m = 0) exists, while in region D, only the two symmetric solutions 
remain. In region B, between the lines T s and T^, the symmetric solution m = is the 
most stable, while in region C, between and T c , the most stable solutions are the 
non-zero ones. The lines T s , and T c meet at a value R* = a/3 - Ik 0.732. 

The nature of the discontinuous transition however calls for a short comment. First 
of all, we note that the evolution equation ( TTTj) is deterministic, so that the system 
should remain trapped forever once it has reached a locally stable stationary solution. 
However, in a more realistic description, the finite size of the system should result in 
a small effective noise, making the system explore the vicinity of the stationary state. 
In this case, the most stable stationary states correspond to the absolute minima of 
gT } R(m) (either m = 0, or two symmetric values m and — m ). 

4. Comparison between mean-field and numerical results 

4-1. Measure of the order parameter 

We now turn to the results of numerical simulations of the agent-based model, and 
compare them with mean-field results. The numerical value for m is obtained as a 
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time average (over 10 4 To) for a simulation usually involving N a = 500 agents and as 
many products as agents (N p = N a ). The memory parameter is set to a = 0.8, a value 
sufficiently large for the continuous time description used in Eq. (|6]) to be approximately 
valid. We have however checked that setting this parameter in the range 0.5 < a < 0.95 
only changes the numerical value of m by less than a few percent for most values of R 
and T. 

We first plot on Fig. H]the steady-state value of the order parameter m as a function 
of temperature T and parameter R. An overall good agreement is observed with the 
mean-field predictions, that are obtained as follows. The fixed points of Eq. (TT7|) are 
solutions of grfii'm) = 0. Besides the trivial solution m = 0, a solution mo > exists 
below the transition line. An analytical expression for m (T, R) is hard to obtain, but 
one can instead easily express T as a function of m and R: 

T(m, R) = K ) l ~ m) (23) 

ln(l + m) — ln(l — m) 

Fixing R, it is then straightforward to plot T as a function of m. One can also invert 
T(R, m) to get R(T, m), as shown on Fig. HI 

4-2. Determination of the transition line 

We also compare on Fig. |6] numerical results for the phase diagram with the mean-field 
predictions. To determine the transition line, we have distinguished two methods for 
the cases of continuous and discontinuous phase transitions. For the discontinuous case 
(R < R*), the transition point can be determined from the jump of the order parameter 
m. In the continuous case, the transition point is more difficult to determine precisely 
from the average value of the order parameter, due to finite size effects. We thus have 
used the standard technique of the Binder cumulant, which characterizes the deviation 
from Gaussianity of the probability distribution of the order parameter [17] . For a 
centered variable m, the Binder cumulant is defined as 

B=l-#4- (*) 

6{m z ) z 

The main interest of the Binder cumulant is that when plotted as a function of the control 
parameter, curves corresponding to different system sizes intersect at the critical point, 
with only small finite size effects. This method allows for a rather precise estimate of 
the transition point. An illustration of the behaviour of the Binder cumulant is given 
on Fig. [51 

Altogether, a rather good agreement is obtained between mean-field and numerical 
results, though some systematic deviations appear, as seen on Fig. El These deviations 
might be due to underevaluated finite size effects (the number of agents used in the 
simulation is N a ~ 500), but it is more likely to result from the simplifying assumptions 
made in the derivation the mean-field equations. Note however that the ranges in 
which continuous and discontinuous transitions are observed numerically are rather well 
reproduced by the mean- field argument. 
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Figure 5. Left: behaviour of the Binder cumulant B as a function of temperature T 
for different sizes N a = 16, 32, 64, 128, 256, 512 (R = 2, simulation time: 5 x 10 5 r ). 
The critical point is obtained as the large N limit of the intersections of successive 
curves. Right: crossing point of two successive curves N a = N/2 and N a = N, as a 
function of N. 




0.3 0.4 0.5 0.6 0.7 0.8 



Q| 1 1 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 4 

R 

Figure 6. Phase diagram (T-R), showing the transition line between the symmetric 
phase and the symmetry-broken phase. Lines are mean-field predictions, and dots 
correspond to the phase transition observed numerically, from the order parameter 
curve for R < R* w 0.732 (discontinuous transition) and from the Binder cumulant 
for R > R* (continuous transition). The mean-field transition line (either continuous 
or discontinuous) is indicated by a full line; it corresponds to T c for R > R* and to T4 
(set of points for which the three minima of the free energy have the same value) for 
R < R* . In the region R < R* , T c is indicated as a dashed line. The inset shows a 
zoom on the region where T c and Td intersect. 



5. Discussion 

In this section, we briefly discuss the relation between the present model and one of 
the most standard models exhibiting a symmetry breaking phase transition, namely the 
Ising model, before drawing some conclusions and perspectives. 
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5.1. Relation to the Ising model 



In the present model, each agent has to choose between two options, namely going 
to the first or the second store. Such a binary choice can be encoded into a 'spin' 
variable s,, equal to 1 if agent i chooses the first shop, and to —1 in the opposite case. 
Though somehow artificial, this definition has the advantage that the order parameter 
m introduced in Sec. 13.31 can be reexpressed as m = N' 1 ^\ Sj. This naturally suggests 
that the present model may exhibit, to some extend, an analogy to the Ising model. 
This analogy can be made more precise, at the mean-field level, in the following way. 
Let us modify the model in such a way that the freshness of products (which drives the 
evolution of the satisfaction of agents) is replaced by the opposite of the energy -E^ing 



(25) 



of the mean-field Ising model, with £r s i ng = — |iV a Jm 2 |j. Then the mean satisfaction a 
introduced in Sec. El evolves according to 

a{t + r ) = aa(t) - (1 - a)AE, 

with 



AE 



m H 

Nn 



m 



1 



-2Jm. 



(26) 



Note that the last equality has been obtained in the large N a limit. We then get the 
evolution equation, in the large memory limit l-a<l, 



da , 

Tt =l{ - 



-a 



2Jm], 



(27) 



which is the analog of Eq. ([6]), with here again 7 = (1 — a)/r . It is easy to check, using 
relation ([7]) to eliminate the variable a, that the time evolution equation for m is given 
by 



dm 
~dt 



2 (1 



m 



-In 

2 



m 



1 + m 



+ Jm 



(2? 



Similarly to Eq. (117)1 . the term between brackets can be written as the opposite of the 
derivative of a free energy f(m), 



T T 1 

f{m) = — (1 + m) ln(l + m) + — (1 — m) ln(l — m) — - Jm 2 



(29) 

2" /v ' 2 K /v '2 y ' 

which is, up to an additive constant Tin 2, nothing but the standard mean-field free 
energy of the Ising model. 

However, beyond mean-field, it is clear that the relation to the Ising model is at 
most an analogy. Indeed, the fact that the interaction between agents is here mediated 
by the ageing of products in each store rather than by an internal energy function makes 
the present model quite different from standard equilibrium models. 

% Strictly speaking, the magnetization m should be computed from the last choice of the N a — 1 other 
agents. 
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5.2. Conclusion and perspectives 

In this paper, we have studied a simple model of competing shops selling fresh products, 
retaining as main ingredients the ageing of products and the satisfaction of customers 
with respect to the freshness of products -considered to be the main factor influencing 
choices made by customers. Under these very simplifying assumptions, we found a 
surprisingly rich dynamics, with both continuous and discontinuous transitions between 
symmetric and symmetry-broken states. The overall agreement between numerical 
simulations and mean-field predictions was found to be satisfactory. 

Several extensions could be interesting to study. For instance, other functional 
forms could be considered for the ageing function h(r) introduced in Eq. (CQ), like 



where 9{x) is the Heaviside function. Such a form would be rather natural for products 
that can be sold only up to an expiration date. We have checked that the form ( 13D1 
preserves the qualitative behaviour of the model. Another natural extension of this 
work would be to include more than two shops. Although we did not study this case 
in detail, a simple argument can already be given here. With a fixed overall number 
of agents N a , adding shops reduces the nominal number of customers for a given shop. 
As a first approximation, a model with more than two shops can be described as a 
two-shop model with an effective number of agents N° s < N a , resulting in an effective 
parameter R eS > R. If one starts from a low value of R corresponding to the symmetric 
region, increasing the number of shops would then eventually lead to a value R in the 
symmetry broken phase. It is thus likely that, for a given total number of agents, there 
is a maximal number of shops that can maintain a good level of activity, which seems 
to be consistent with conventional wisdom. 

Further possible extensions may include taking into account explicitly the prices 
of products, as well as the costs associated to the shops (e.g., rent or salaries). In this 
context, one can imagine to include, in a simplified way, a dynamics of the shopkeeper 
who could modify the prices, or diversify the products, as a function of the dynamical 
state of the store. One might guess that such a dynamics could contribute to stabilize 
the symmetric states. 
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